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Abstract 

We present a theoretical reappraisal of the branching ratios and CP asymmetries for the decays 
B — > X q £ + £~, with q = d,s, taking into account current theoretical uncertainties in the description 
of the inclusive decay amplitudes from the long-distance contributions, an improved treatment of 
the renormalization scale dependence, and other parametric dependencies. Concentrating on the 
partial branching ratios AB(B — > X q £ + £~), integrated over the invariant dilepton mass region 
1 GeV 2 < s < 6 GeV 2 , we calculate theoretical precision on the charge-conjugate averaged partial 
branching ratios (AB q ) = (AB(B -► X q l+£-) + AB{B -> X q l + l~))/2, CP asymmetries in partial 
decay rates (a CP ) q = (AB(B -> X q £+£-) - AB(B -> X q £+£-))/(2(AB q )), and the ratio of the 
branching ratios A1Z = (ABd) / (AB S ) . For the central values of the CKM parameters, we find 
(AB S ) = (2.22lg;») x 1(T 6 , (AB d ) = (9.6lt};f 7 ) x 1(T 8 , (a CP ) s = -(0.19±°;^)%, (a CP ) d = 
(4.40±|;|^)%, and ATI = (4.32±0.03)%. The dependence of (AB d ) and ATI on the CKM parameters 
is worked out and the resulting constraints on the unitarity triangle from an eventual measurement 
of A7Z are illustrated. 
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1 Introduction 



With the advent of new and upgraded experimental facilities in the next year(s), flavour physics 
involving B decays will come under minute experimental and theoretical scrutiny. The overriding 
interest in these experiments is in measuring CP-violating asymmetries in partial -B-decay rates, 
which will allow to quantitatively test the Kobayashi-Maskawa |Qj paradigm of CP violation. In 
addition, the large number of B hadrons anticipated to be produced at these facilities (estimated 
to be O(10 8 ) - O(10 12 )) will allow to measure a number of flavour-changing-neutral-current (FCNC) 
processes involving the transitions b — > sX and b — ► dX, with X = 7, g, £ + £~, vv, and B° - B° 
mixings. In the context of the Standard Model (SM), FCNC decays and mixings measure the Cabibbo- 
Kobayashi-Maskawa (CKM) Q matrix elements, in particular Vtd, Vt s and Vtb- These quantities can, 
in principle, also be measured directly in top quark decays t — > qiW + , with = d, s, b. A comparison 
of these matrix elements in the FCNC processes and direct measurements in t decays would provide 
one of the best strategies to search for new physics in B decays. So far, only Vtb has been directly 
measured at Fermilab, yielding \Vtb\ = 0.99 ± 0.15 ||]. 

Present knowledge of Vtd owes itself to the measurements of AMd, the mass difference in the B° - 
B° complex. With the current world average AM d = 0.471±0.016 (ps)" 1 , the error on Vtd is dominated 
by theoretical uncertainty on the hadronic matrix element fs d y/B^, for which present Lattice-QCD 
estimates are f Bd ^/B^ = 215 ± 35 MeV |], yielding 0.0065 < \V td V* b \ < 0.010. We also mention 
that a single event for the charged kaon decay mode K + — > ix^vv reported by the Brookhaven E787 
experiment, yielding B{K+ tt+vv) = (4.2+§;^) x 10~ 10 , allows one to infer 0.006 < \V td V^\ < 0.06 
The branching ratio for the decay B — > X s j has led to a determination of the matrix element 
Vt s fj) yielding |Vf S V^|| = 0.0035 db 0.004, with the error dominated by the experimental error on the 
branching ratio B{B — > X s + 7) [||,|7]]. These numbers can be taken as the measurements of \Vt d \ and 
I Vf S I by assuming the value Vtb — 1 from the CKM unitarity, which holds to a very high accuracy ||]. 

In this paper, we pursue the idea of measuring the FCNC semileptonic decays B — > X s £ + l~ and 
B — > Xd£ + £~, below the J/ip- and above the p, w-resonance regions in the dilepton invariant mass, to 
determine \Vt s \ and \Vtd\, respectively, and the ratio |Vtd/Vt s | from the ratio of the branching ratios. In 
this context, these decays and the related ones, B — * X s vv and B — * Xdyv, were discussed some time 
ago |^]. The decays B — ► (X s , Xd)vv are practically free of long-distance complications |[0| and the 
renormalization-scale dependence of the decay rates has also been brought under control 
these decays are theoretically remarkably clean but, unfortunately, they are difficult to measure in 
T(4if?) decays and out of question in hadronic collisions. Using the missing energy technique and LEP 
I data, the ALEPH collaboration has searched for the decays B — » X s vv setting an upper bound 
B(B -> X s vv) < 7.7 x 10~ 4 (at 90% C.L.) Q, which is a factor 20 away from the SM expectations 
[pd|]. While the discovery of these decays looks formidable elsewhere, a high luminosity Z°-factory - 
being discussed in conjunction with an e + e~ linear collider looks like having the best chance of 
measuring them. This possibility deserves a dedicated study. 
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The possibility of determining |Vfd/Vt s | from the ratio of the invariant mass decay distributions 
= ^[B — > Xd£ + £~]/^Jj[B — ► X s £ + £~\ away from the resonances was revisited by Kim, Morozumi 
and Sanda ||14j]. These authors included the effects of the leading order power corrections (in 1/m 2 ) in 
the short-distance part of the dilepton invariant mass distribution and the long-distance contributions 
from the cc-resonances, calculated in Ref. [15|. (For earlier- vintage derivations without the power 
corrections, see [16,17|.) We reanalyze the decays B — > X s £ + £~ and B — > and the ratio of 

the branching ratios AT?. = / ds|f[.B — > X d £ + £~}/ J ds^[B — > X 8 £ + £~], integrated over a kinematic 
range < s < q^ax^ designed to minimize the resonant contribution. Our theoretical treatment 



differs from that of Ref. [14] in a number of ways, summarized below. 



The dilepton invariant mass distributions in B — > (X s , X d )£ + £~ can be calculated in the con- 
text of the heavy quark effective theory (HQET) as a power expansion in regions far from 
the resonances, thresholds and end-points [jn^JlOfl . Away from the J/tp, tp', ...-resonances, the 
1/m^-expansion provides, in principle, a viable description of the non-perturbative contribu- 



tions arising from the cc-loop [10|. The contribution of the light quark (/(/-loops, which is not 
CKM-suppressed in the decay B — ► X d £ + £~ , can likewise be calculated by doing an expansion of 
the decay amplitudes in A.Q CD /q 2 in regions of the dilepton squared mass satisfying q 2 3> Aq CD . 
Thus, the HQET framework provides an evaluation of the invariant dilepton mass spectrum in 
these processes with the present precision limited to the leading power corrections in 1/ m 2 , 1 /m 2 
and AQ CD /q 2 . We present HQET-based calculations of the decay rates, CP asymmetries and 
the ratio A1Z. 

Away from the resonances and the end-points, the power corrections in 1/m 2 calculated in 
HQET and in explicit wave function models, such as the Fermi motion (FM) model |l8ft , yield 



very similar invariant dilepton mass [15] and hadron energy distributions |ll| in the decays 
B — > X q £ + £~ . However, it is known that there are marked differences in estimates of the non- 
perturbative cc-contribution, obtained by using the 1/m^-corrections in the HQET approach 



and alternative methods based on the Breit-Wigner-shaped resonant amplitudes p0|,21|. Data 
may eventually provide a discrimination against some of these approaches, but currently at least 
four different variations on this theme exist in the literature [10, 15,^2, 23|. This LD-uncertainty 



therefore compromises theoretical precision on decay rates and has to be taken into account. We 
calculate the theoretical uncertainties on the branching ratios for the decays B — > (Xd, X s )£ + £~ , 
CP asymmetries and the ratio A1Z, showing numerically their impact on the determination of 
\Vts\, \Vtd\ an d the CKM-Wolfenstein parameters p and r/ |24| from an eventual measurement of 
these decays. 

We reanalyze the renormalization scale dependence in the branching ratios for the decays B — * 
X s £ + £~ and B — > X d £ + £~ , using the method employed by Kagan and Neubert in the radiative 
decay B — ► X s +j [p5]| . This approach avoids accidental cancellations among the individual scale- 
dependent contributions but gives a larger scale (^-dependence of the branching ratios than the 



method of evaluating the same in the total branching ratio [O] . The former is probably a more 
realistic estimate of the neglected higher order corrections. 

We find that the partial branching ratio in the SM is uncertain by typically ±13% (±15%) for the 
decay B — ► X s £ + £~ {B — ► Xd£ + £~), but the ratio A1Z is remarkably stable with typical error less 
than several percent. Hence, A1Z is well-suited to determine the ratio |V^/Vf S |. However, the scale- 
dependence of the CP asymmetries in B — > (X s , Xd)£ + £~ is found to be huge, reflecting the (present) 
leading logarithmic theoretical accuracy of the CP-odd parts of the amplitudes. Without the power 
corrections and fixing the scale to fi = m&, the CP asymmetries in question have been studied earlier 
in Ref. [26]. We point out that these estimates are uncertain by almost ±100% due to the sensitive 



scale-dependence and their stabilization requires next-to-leading order corrections. In the case of the 
CP-even parts, we recall that the inclusion of the explicit 0(a s ) corrections in the matrix elements 
has reduced the scale dependence of the decay rates considerably [ 27,28)]. 



This paper is organized as follows: In section 2, we briefly review the derivation of the matrix 
elements and dilepton invariant mass distributions for the decays B — ► (X s ,Xd)£ + £~ including long- 
distance contributions in the four approaches: (i) AMM fT^JH, (ii) KS f22j, (iii) LSW [§3| and (iv) 



HQET [ 10 1 . The partially integrated branching ratios and CP asymmetries are presented in section 
3 where we also specify our input parameters. We show the scale dependence of the branching ratios 
AB(B — » X s £ + £~) and AB(B — > Xd£ + £~) in the AMM approach and the contributions arising from 
the individual Wilson coefficients. We also present a comparative numerical study of the quantities 
(AB S ), (ABd), (acp)s and (acp)d hi the four mentioned approaches. Uncertainties arising from the 

(6) 

other parameters (m^, m< and h-q CD ) are worked out numerically. With this we calculate the overall 
theoretical errors in these quantities and the ratio A1Z and their impact on the determination of the 
CKM parameters. Finally, section 4 contains a brief comparison of the theoretical precision of | Vtd/Vt s \ 
in the decays B — > (X s ,Xd)£ + £~ with that of other methods proposed in the literature to determine 
the same ratio. 

2 B — > (JQ, X s )£ + £~ Decays in the Effective Hamiltonian Approach 

We work in the effective Hamiltonian approach, which is based on integrating out the heavy degrees of 
freedom (t, W^, Z°), in the SM. The resulting effective Hamiltonian for the decays B — > (Xd, X s )£ + £~ , 
TL e ff(b — » q£ + £~), can be expressed as follows: 

n eff (b - q£ + n = -^-VtlVt^dOi + ^=-v: g v ub [ci(Oi< u > - Oi) + c 2 (o 2 w - o 2 )] , (l) 

where Vij are the CKM matrix elements. The C{ are the Wilson coefficients, which depend, in general, 
on the renormalization scale /U, except for Cio, and can be seen in leading logarithmic approximation 



in 1 27 ] . The operators are defined as follows: 

Oi = {qL a l^ La )(c L pYc L p) , 
3 



2 = (qLa7 f ib Lf 3)(c Lf3 j fl c Ija ) , 

C 6 = (qLal^Lp) (.^Rfi'f^Ba) > 

q f =u,d,s,c,b 

e 2 

° 9 = TTT^qal^Lb^^ , 



e 2 
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^fcT^WT^ , (2) 



where L and i? denote chiral projections, £(-R) = 1/2(1 ^75). Here, unitarity of the CKM matrix has 
been used in writing the flavour structure of a generic FCNC 6 — > g transition amplitude T'W in the 
form 



r (,) = £ A^ = AS 9) (^-T c ) + AW(T u -r c ), (3) 



i=it,c,i 



where A^ = V*Vib and q = d,s. For the 6 — ► s transitions, the second term in Eq. (^) can be safely 
neglected as Atf <C A( S \ However, for the b — > d transitions, the CKM factors Au and A^ are of the 
same order and hence all terms in Eq. (J3j) must be kept. The operator basis given in Eq. (1) has been 
written in accordance with Eq. (|3|) and includes the Four-Fermi operators containing a uu pair, 

Oi {u) = {qL a ln b L*){uLf3l^U LI3 ) , 

2 {u) = (q La ^b Lf3 )(u L ^u La ) . (4) 



The matrix element for the decays b — > ((7 = d, s) can be written as 

M(b^q£ + r) = ^V t * q V tb [(C^ - C 10 ) (q^Lb) L£) + (C^ + C 10 ) (q^Lb) (I^ R£) 

(5) 



V2tt 

- 2Cf (q i <T M „ ^ (m g L + m 6 i?) ^ (£y £) 

Here g v = p+ + pi denotes the Four-momentum of the invariant dilepton system, where p± are 
the corresponding momenta of the £^\ s = q 2 is the invariant dilepton mass squared. The effective 
coefficients of Og are given by 

C$(s) = C 9 r)(s)+Y«(s). (6) 
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The functions r/(s) and Y q (s) represent the 0(a s ) correction [^] and the (perturbative) one loop 
matrix element of the Four-Fermi operators |27],^], respectively. We have in the (naive dimensional 
regularization) NDR-scheme, which we use throughout our work, 



Y«(s) 



g(m c ,s) (3d + C 2 + 3C 3 + C 4 + 3C* 5 + C 6 ) 

s) (4 C 3 + 4 C 4 + 3 C 5 + C 6 ) - i 5 (0, a) (C 3 + 3 C 4 ) 
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(3C 3 + C 4 + 3^ + 00 



(3Ci + C 2 )( 5 (0,s)- 5 (m c ,s)) , 



(7) 



where we have introduced the dimensionless variable s = q 2 /m 2 and m c = m c /mb. The functions rj(s) 
and g(z,s) can be seen elsewhere [ p7| , p0| l. Note that the renormalization scheme-dependence of the 
function Y 9 (s) cancels with the corresponding one in Cg. The effective coefficient of the bs^f vertex is 
given by Cf ff = C 7 - C 5 /3 - C 6 |5|. 

The dilepton invariant mass spectrum including power corrections in the HQET approach in B — ► 
X q £ + l~ decays can be written as: 



dB 

ds 



dB° dB 1 ' 1 



ds ds 

where the first term corresponds to the parton model 



+ 



dB 1 /* 2 



ds 



(8) 



pq| , the second term accounts for the 0(1/ ml) 
power corrections [15|, and the last term accounts for the non-perturbative interaction of a virtual uu- 
and cc-quark loop with soft gluons. The explicit expression for dB x l q2 jds for m q = can be deduced 



from the literature [1C] 



ds 



-£ C 2 A 2 p(l 



(9) 



x Re 



G 



eff* (! + 6 « 



+ C^°>(s)(2 + s) 



F(s,m c ) 



\(«) 



F(s,m u ) F(s,m c ) 



III,: 



+ [(3Ci + C 2 )(g(0, s) - g(m c , s))}* (2 + s) 



,A^ 2 F(s,m u ) 
•(g) 1 l " 



F(s,m c 



Ai 9) F(s,m c 



mt 



mt 



A 



mi 



The branching ratio for B — > X„£ + £ is expressed in terms of the measured semileptonic branching 



ratio B s i for the decays B 



X c £i/£. This fixes the normalization 

3a 2 \V t *V tb \ 2 



B = B S 



1 



(10) 



iGvr 2 \V cb \ 2 f(m c )K(m c ) ' 

where f(m c ),K(m c ) can be seen, for example, in fl5|| . The function F(s,m) = F(r) with r = s/(4m 2 ) 
is given in flOf . In the region r> 1, F(s,m u )/m 2 l oc 1/s. The condition r 1 is well satisfied, for 
example, for q 2 > 1.0 GeV 2 (for which r > 25). In this region, the operator product expansion (OPE) 
is not in 'l/m 2 ' but in I& CD /q 2 . Hence, there is a sufficiently large region in q 2 where the OPE holds 
in 1/m 2 , 1/m 2 and XQ CD /q 2 . Note also that for the terms proportional to the power corrections, we 



use 



Cgf {0) *(s) which equals C$*(s) with r?( 



9q 



1. 

5 



In B — > X g ^ + £~decays cc-resonances are present via B — > Ag + (J/ip, ifi', ...) — > X q £ + £~. Their 
implementation and the corresponding uncertainties in the -B — ► X s £ + ^~ case have been discussed 
recently by us |20|. There are at least four different Ansatze advocated in the literature in this 
context, summarized below. 

• The HQET-based approach [fT0|], where the non-perturbative cc-contribution away from the 
(J/ip^tp 1 , ...) -resonances is implemented by the 1/mj! terms in the expression for dB 1 ^ 2 /ds. 

• One could add the resonant cc-contribution, parametrized using a Breit-Wigner shape with 
the normalizations fixed by data, to the complete perturbative contribution resulting from the 

-loop. This scheme has been used in a number of papers p~7 Jl5 , 14 2C ] . 



cc- 



The effective coefficients including the cc-resonances are defined as 

Cff (§) = C 9V (s) + Y*(s) + Y res i(s) , (11) 
where Y q (s) has been given earlier and Y res q (s) in this scheme is defined as: 

x £ , (12, 



with C<°> s 3d + C 2 + 3(7 3 + C'i + 3C 5 + C„. In what follows we shall neglect the part 



v t * q v tb 

in Eq. (|l^) in our numerical analysis, since the particular combination of the Wilson coefficients 
appearing in this term is strongly suppressed compared to C^°K Further, since data only determines 
the product kC^ = 0.875 we keep this fixed. For ease of writing, we call this approach the AMM 
approach |L7]]. 

The remaining two approaches are the following: 

• The LSW-approach p3[| : Here, for the non-resonant cc-contribution, only the constant term in 
g(rh c , s) is kept. Calling it g{rh c , s), it is given by g(rh c , s) = — | ln(mf,//u) — | lnm c + ^. The 
resonant cc part is essentially as given in Eq. (|l^). 

• The KS-approach [p^j , in which the function Cg^(s) is parametrized using a dispersion approach. 
For details and further discussions of this approach, we refer to [p2|,p0|. 



In B — > Xdl + £~ decays, in addition to the cc bound states, also the uu bound states have to be 
included in the decay amplitudes. We have calculated the dilepton invariant mass distribution, using 
the Breit-Wigner shape for the resonances, as discussed earlier, and taking the widths and partial 
leptonic widths from the Particle Data Group However, numerically the uu-resonant part is less 
important, as the leptonic branching ratios B(V° — > e + e~) and B(V° — > fj, + fi~) for the dominant 
resonances V° = p°,co are small 0. Moreover, their effect is reduced by imposing a cut on the 



dilepton invariant mass, say q 2 > lGeV 2 , which we have explicitly checked. Higher states like p',ui' 
have larger widths and are thus expected to play minor roles due to their smaller branching ratios in 
dilepton pairs. 

In the three approaches discussed above (AMM,LSW,KS) we include the l/m 2 -corrections, calcu- 
lated in the phenomenological Fermi motion model (FM) , which implements such effects in terms 
of the £>-meson wave function effects. The implementation of the FM model in B — > X s i + l~ decays in 
the dilepton invariant mass distribution can be seen in J15[| , which we also adopt here for the calcula- 
tions of the distributions in B — > X^l^. We note that the branching ratios in the HQET-based 1/m 2 
approach and the FM-model are very close to each other for identical values of the input parameters. 

3 Branching Ratios and CP Asymmetries in B — > X q £ + £~ 

3.1 Numerical input and definitions of the partial branching ratios and CP asym- 
metries 

We now specify how we determine theoretical uncertainties in the branching ratios, the ratio A1Z, 
and CP asymmetries in the decays B — > (X s , X^l+i - . The dispersion in the values of the obsevables 



my/ 


80.41 GeV 


m z 


91.1867 GeV 


sin 2 9w 


0.2255 


m s 


0.2 GeV 


m d 


0.01 GeV 


m b 


4.8 ± 0.2 GeV 


m t 


173.8 ± 5.0 GeV 


a(5) 
1V QCD 


0.220lffi GeV 


a- 1 


129 


a s {mz) 


0.119 ±0.0058 


B s i 


(10.4 ± 0.4) % 



Table 1: Default values of the input parameters and the ±1 a errors on the sensitive parameters used 
in our numerical calculations. 

(5) 

due to the errors in the input parameters mb,fi,m t , a s (mz) (equivalently Kq CD ), and B s i, given in 
Table |l|, is calculated by varying one parameter at a time. To estimate the uncertainty from the 6-quark 
mass in the FM model, we explore the parameter space of this model with three sets of parameters: 
(pF,m q ) = (520,280), (450,0), (245,0) in (MeV,MeV), which correspond to an effective 6-quark mass 
of mf^ = 4.6, 4.8, 5.0 GeV, respectively. We set m c = m^(m^) — 3.4 GeV in both the FM-model 
and HQET analysis. Comparison with the HQET prediction |l5| is worked out for Ai = —0.20 GeV 2 
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and A2 = 0.12 GeV 2 , as the dependence of the branching ratios on these parameters is small. The 
individual errors are then added in quadrature to get the final cumulative error. 
We proceed by defining the partly integrated branching ratios (q = s,d): 

2 dB(B -» X q £+£-) 



AB n 



q 2 

imax 



2 

min 



dq 1 



dq 2 



(13) 



together with AB q , for the CP-conjugate decays B — > X q £ + £ , and the branching ratio averaged over 
the charge-conjugated states: 



(AB q 



AB a + AB a 



The CP asymmetry in the partial rates for B —* X q £ + £ is defined as: 



AB q - AB q 
AB q + AB q 



(14) 



(15) 



We further decompose the partial branching ratios AB q in terms of the CKM factors 

AB q = (|Ah 2 D?> + lA^I 2 ^) + Re(X^X^)DM + I m {\^X^)D^)/\V^ , 

from which the CP conjugated branching ratio AB q can be obtained by substituting X^\ 
Hence, the charge-conjugate averaged branching ratio (AB q ) is obtained from AB q by dropping the 
Im(x[ q ^* Xu^) term. The CP asymmetry is given by the expression: 



(16) 
v(g)* 

Vi ■ 



(a CP ) q = Im(x{ q) *X^)Dl q) /(\V cb \ 2 (AB q )) . 



(17) 



The functions , j = t,u,r,i depend on the input parameters, which we have specified in Table 1, 
and on the interval in q 2 , specified by q^n and q^ax- We shall work always above the (p, lo)- and 
below the J/V'-resonances in the so-called low-g 2 region with q^^ and q^ax taken as 



7min 



1.0 GeV 2 < q 2 < 6.0 GeV 2 



7max 



(18) 



W) 



-AX 
AX 2 (1 



A 2 



We use the Wolfenstein representation of the CKM matrix \24\ with A = 0.819 and A = 0.2196 
fixed, as the errors on these quantities are small ||. The other two parameters (p, rj) are implicitly 
the subject of the present work. Defining p = p(l — 4r) and fj = rj(l — 4^), we have up to terms of 
order A 6 gl): 

AW = AX A (p - i-q) , AP' - ---1A- 1 

x^ = Ax\- P -in) , x 

and V cb = AX 2 . It follows that ^ 2 = A 2 (l + A 2 (l - 2p))((l - p) 2 + f) + 0(A 6 ). Global fits of the 
CKM parameters have been performed in a number of papers |p2| , p3| ) p4]| , with very similar (though 
not identical) results. For illustration, we shall use the results of the CKM fits from Ref. |32|], yielding: 



2 

p + iff) , 



+ X 2 (p-irj) 



(19) 
(20) 



0.155 



-0.115 
-0.105 ' 



V 



0.383 



-0.063 
-0.060 • 



(21) 
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-I— - I 



6 

/j.[GeV] 



10 



6 

^[GeV] 
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Figure 1: Renormalization scale (fj,)- dependence of the individual terms and the partly integrated 
branching ratios AB S for the decay B — » X s £ + £~ (a) and ABd for B —* Xdl + £~ (b), calculated in the 

AMM-approach. The solid, dotted, dashed, long-short dashed curves correspond to the contributions 

eff eff eff eff 

proportional to the effective Wilson coefficients \C 7 | 2 , |Cio| 2 , |C 9 | 2 and Re(C 7 C 9 ), respectively. 

The resulting /i uncertainty in the branching ratio, obtained by adding the weighted errors in quadra- 
ture, is indicated by the shaded area. 

3.2 Parametric dependence of the branching ratios and CP asymmetries 

We study the scale (/^-dependence of the branching ratios along the lines followed in |2[| in the 
B — > Xsj case. Thus, instead of varying the scale \x between and Im^ in the full expression for the 
respective branching ratios (the naive method), the scale-dependence of the individual terms involving 
different Wilson coefficient combinations is calculated independently and the resulting errors are added 
in quadrature. It is a conservative approach and avoids the possibility of accidental cancellations of 
the scale-dependence in the various terms, which takes place in the SM in both the B — > X s j case 
] and in B — > X q £ + l~, as shown here. For the branching ratio in B — > X g ^ + £~decays the relevant 
coefficients are: |Cio| 2 , |Cg^| 2 , Re(Cf^Cg) and |C*7^| 2 . Of these, C%o does not renormalize, however, 
there is a residual dependence on \x from the normalization for which inclusive semileptonic branching 
ratio is used, bringing in an extra a s (/^-dependence. 

The scale-dependence of the individual contributions from the specified Wilson coefficients to the 
branching ratios AB S and ABd and the branching ratios themselves, are shown in Fig. 1(a) and 1(b), 
respectively. We find for the scale dependence of AB S an uncertainty (+9.0, —7.3)%, measured from the 
reference value [i = m,b- This is to be compared with the corresponding uncertainties (+4.1,-1.3)% 
calculated in the naive approach. The estimated /i-dependent uncertainty in ABd 1S found to be 
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(+7.7, —7.6)%, compared to 2% in the naive approach. 
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ID n3 


4 90 


m t = 178.2GeV 


2.36 


-0.18 


10.18 


4.18 


m t = 168.2GeV 


2.10 


-0.20 


9.06 


4.63 


^qcd = 0-298GeV 


2.20 


-0.16 


9.52 


3.74 


^qcd = 0.157GeV 


2.24 


-0.22 


9.70 


5.03 



Table 2: Values of the charge- conjugate averaged partial branching ratios (AB S ) and (ABd) and the CP 
asymmetries (acp) s and (acp)d, in the four LD- approaches AMM JT^/ , KS LSW $2^J and HQET 
J7^/ , discussed in the text. In the top part of the table (above the horizontal line), the parameters are 
fixed to their central values given in Table 1 and Eq. In the lower part of the table, the parametric 

(5) 

dependence of the observables on m^, mt and Kq CD , calculated using the AMM-approach, is listed. 

The dependence of the charge-conjugate averaged branching ratios (AB S ) and (ABd), and the CP 
asymmetries (acp) s an d (acp)d on the four schemes concerning the cc-contribution is shown in the 
upper part of Table 2. For all these entries, we have fixed the parameters to their central values given 
in Table 1 and Eq. (|2l|). The dependence of these observables on m^, mt and Aq^ cd , obtained in the 
AMM-scheme by varying only one parameter at a time, is shown in the lower part of Table 2. For 
the central values of p and r/, the partial branching ratios are found to vary in the four approaches in 
the range: 2.05 x 10~ 6 < (AB S ) < 2.31 x 10~ 6 and 8.83 x 10~ 8 < {ABd) < 9.98 x 10~ 8 . For the same 
values of p and r\ but taking into account in addition the rest of the parametric uncertainties in Table 
2, B s i, and the scale-dependence from Fig. 1(a) and 1(b), we find: 

(AB S ) = (2.22±gjg) x 10- 6 , 

(AB d ) = (9M±\%) x 10~ 8 . (22) 

Thus, apart from the CKM-parametric dependence, we estimate ±13% uncertainty on {AB S ) and 
somewhat larger, ±15%, on (ABd). These errors are significantly larger than what one comes across 
in the literature. The present experimental bound is B(B — » X s l + £~) < 4.2 x 10~ 5 (at 90% C.L.) 
[]35| 1. We are not aware of a corresponding bound on B(B — > Xdl + £~). 

The branching ratio (ABd), calculated in HQET, is shown in Fig. |2| as a function of the CKM 
parameter p for three fixed values of rj, which correspond to the central value and the 95% C.L. 



bounds given in Eq. (21). The other input parameters have been fixed to their central values given 
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in Table |]. In the allowed CKM parameter space, this partial branching ratio varies by a factor 3. 
As the theoretical error from the rest of the parameters is estimated to be ±15%, the measurement 
of (ABd) should allow to determine p and rj. The ratio A1Z = (ABd) / {AB S ) has lot less theoretical 
error, as shown below. 



The CP asymmetry, (acp) s defined in eq. (15) in the b —* s case in the SM is small. Hence its 



measurement can be used to search for new sources of CP violation in the b — * s£ + £~ transition. 
Numerically, the CP asymmetries are more uncertain reflecting in particular the scale-dependence of 
the functions . A qualitatively similar behaviour has also been noted for the CP asymmetries 



in the radiative decays B — > X s + 7 and B — ► Xd + 7 in |37y . However, the scale-dependence of 
the CP asymmetries is more marked in the decays B — ► (X s , Xd)f r t~ due to cancellations in two 
different products of the Wilson coefficients entering in . (Specifically, between Cf^Im(Cg^\ u ) 
and Im(Cgg\ u Cgg\t), with Cgjf\ x denoting the part in which is proportional to the CKM factor 
This can be seen in Fig. ^, where we show the //-dependence of the two mentioned contributions 
in d\ , and the function itself calculated in the naive and independent approaches. The function 
is very similar and hence not shown. The //-dependence of in the naive approach, shown 
by the long-short dashed curve, is very marked and it gets further accentuated in the independent 
approach, shown by the two dashed curves. For the central values of the CKM parameters and 
estimating the //-dependence in the independent approach, we find: 

(a CP ) s = -(0.191°;^)%, 

(a CP )d = (4.40^)%. (23) 

The corresponding numbers in the naive scale-dependent method are: (acp) s = — (0-19^ ' 13 )%, an d 
(acp)d = (4.40i§;23)%. ^ n either case, Fig. ^ underscores the importance of calculating the next-to- 
leading order effects in (acp) q - 



3.3 Extraction of 



YtA 

Vts 



For a precise determination of jj^4 (equivalently the CKM parameters), we calculate the ratio: 

A*^. (24, 
In terms of the CKM parameters and the functions d\ s ^ and Dj with j = t, u, r, defined earlier: 

2 , ~2\ n (d) , ( -2 , _2s n (d) , ( - ( , _ -s _ - 2 \ n {d) 



Jl-pf + f)D\ a > + {f + ?)D™ + {p{\ -p)- fj ) D 



ATI = A 2 — r > ' •' ' yr ' •' '— r ' ■' '~ r , (25) 

(1 - A2(l - 2p))D[ s) 

(s) (s) 

where we have neglected terms proportional to Xu /X\ . A simpler form for A7Z follows, if one notes 
that the functions and are equal for all practical purposes (see Table 3). Hence, setting 

D (d) _ D w one has 



(l-A2(l-2p) 



(l-p) 2 + r/ 2 )£» t ^ (l-p) 2 + r/ 2 ) £,(') 
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(26) 



-0.1 0.1 0.2 0.3 0.4 

P 

Figure 2: The charge-conjugate averaged partial branching ratio {ABd) in the HQET-approach for 
the decay B — ► Xd£ + l~ as a function of the CKM parameter p for three values of rj; solid curve 
(n = 0.383), dotted curve {rj = 0.5), dashed curve (rj = 0.27). 



The overall CKM factor is just the ratio \Vtd\ /\Vts\ ■ Note that the first (and dominant) term is 
independent of the dynamical details. The ratio Dy® /d\ s ' is found to be numerically small (but 
model dependent, varying between 1.03 x 10 -2 for the KS-approach and 2.16 x 10 -2 for the LSW 
approach). The ratio /D^ is, in general, larger and it depends more sensitively on the estimate 
of the long-distance cc-contribution, varying between +0.14 (for the LSW-approach) and —0.12 (in 
HQET). However, the multiplicative CKM factor accompanying this term in Eq. (|26| ) being small 
comes to rescue. For example, for p = 0.151 and f\ = 0.374, this factor is only —0.012. Hence, for 
these values, we find A1Z = (4.32 ± 0.03)%. For other values of the CKM parameters, the uncertainty 
is larger and we quantify it later. The ratio A1Z as a function of p is shown in Fig. [4] for the HQET- 
method. The three curves correspond to rj = 0.5 (dotted curve), rj = 0.383 (solid curve), and r/ = 0.27 
(dashed curve). 

We now evaluate the theoretical precision in the determination of 
ment of A1Z. The other uncertainties being insignificant, there are 



lid 



t— from an eventual measure- 

Vts 

basically two sources of errors: 
(i) a small residual scale-dependence, and (ii) the LD-scheme-dependent uncertainty, which depends 
on the parameters p and r\. In Fig. || we show the constraints on p and r\ from an assumed value 
of A7£ with the LD-effects calculated in the AMM-approach. For each value of A1Z, the practically 
overlapping curves represent the effect of varying p in the range < p < 2m&. Numerically, the 

net p uncertainty on the ratio A1Z is found to be ±0.6%. The effect of the errors of mt, a s (mz) and 
the 6-quark mass are smaller and not shown. 

The potentially largest uncertainty in A1Z, due to the LD-effects, is shown in Fig. |(| where we have 
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-6 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

2 4 6 8 10 

fj,[GeV] 

Figure 3: Renormalization scale (/i) -dependence of the individual contributions and the function 
D^, calculated in the AMM- approach. The solid and dotted curves correspond to the contributions 
proportional to the effective Wilson coefficients Cf^Im(Cg^\ u ) and Im(Cg^\ u Cg^*\t) , respectively. 
The naive p dependence is shown by the long-short dashed curve. The resulting p uncertainty in the 
independent approach is bounded by the dashed lines. 

plotted the constraints on p and r\ from assumed values of A1Z. The four curves shown correspond 
to the LD-schemes: AMM, KS, HQET and LSW. As remarked earlier, the LD-related uncertainty is 
vanishingly small for the central values of p and rj, i.e. at or close to the apex of the drawn triangle. 
However, for other points in the (p, ?7)-plane, the uncertainty is perceptible but still small, except for 
regions of the (p, r^)-plane which are already ruled out from the existing CKM fits. 



4 Theoretical Precision on |V^/V^ S | from B Decays 

The ratio A1Z should be measurable at the Tevatron, the later phase of the B-factories, and certainly 
at the LHC. The merit of ATZ lies in the theoretical precision on | Vtd/Vt s \ (or on the unitarity triangle) 
which we have estimated here and found to be quite competitive with other proposals in the market, 
some of which are reviewed below. 

The B°-B° mixing ratio AM S / AM d can be expressed as follows: 

AM, _ M Bs (f 2 Bs B Bs ) V t 2 

AM d M Bd (f Bd B Bd yV td l • 1 ' 

The achievable accuracy on Vtd/Vt s depends, apart from the experimental measurement error, on 
the knowledge of the ratio of the hadronic matrix elements £ = f Bd \JB Bd j f Bs ^JB Ba , for which the 
current Lattice estimate is £ = 1.14 ± 0.06 ± 0.03 ± 0.10 |J. The errors reflect, respectively, the actual 
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^ d, [io- 8 ] 


^ (s) [io- 6 ] 


A (d) [l0" 8 ] 






AMM 


2.31 


3.75 


2.30 


20.96 


-2.34 


-2.34 


KS 


2.12 


2.18 


2.11 


1.42 


-2.00 


-2.05 


LSW 


2.40 


5.16 


2.39 


32.59 


-2.50 


-2.43 


HQET 


2.14 


2.88 


2.13 


-24.89 


-1.99 


-1.94 


m b = 4.6GeV 


2.24 


4.48 


2.22 


26.83 


-2.31 


-2.26 


m b = 5.0GeV 


2.41 


3.47 


2.40 


18.86 


-2.39 


-2.31 


m t = 178.2GeV 


2.45 


3.75 


2.44 


21.89 


-2.36 


-2.35 


m t = 168.2GeV 


2.18 


3.75 


2.17 


21.61 


-2.33 


-2.33 


A® CD = 0.298GeV 


2.29 


3.39 


2.28 


20.71 


-1.97 


-1.95 


AgL = 0.157GeV 


2.33 


4.15 


2.32 


21.35 


-2.70 


-2.73 



Table 3: Values of the functions Dj, j = u,t,r,i and D^\d^ s ' defined in eq. (25) and (17) in the 
four schemes discussed in the text for the central values of the input parameters. The entries below 
the horizontal line correspond to using the AMM scheme, and varying the input parameters, one each 
at a time, fixing the rest to their central values. 



calculational error of this ratio in the quenched approximation, estimated effects of unquenching, and 
from chiral loops. Thus, the present theoretical error on this quantity is of O(10%) and it remains a 
theoretical challenge to improve this significantly. However, the measurement of AM S , for which the 



present experimental lower bound is 12.4 ps _1 (at 95% C.L.)[33], may turn out to provide the first 
measurement of Vtd/Vts, as the central value of AM S in the SM is around 14 ps _1 [ 32 , 33 , 34 1 , which is 
not too far from the present limit. 

Theoretical precision on ATZ is comparable to the one on the corresponding ratio of the branching 
ratios involving the CKM-suppressed decay B — > Xd + 7 and the CKM-allowed decay B — > X s + 7 
p6|,|3"7|] . Defining the ratio of the branching ratios as (implied are charge-conjugate averages) 

(B(B -^X d + 7)) 



R{d 1 )/s 1 ) = 



(28) 



(B(B^X s + 7 ))' 

the ratio R(d r y)/s , j) gives a constraint on the CKM matrix elements which is very similar to the one 
given by ATZ (compare Eq. (26) in Ref. [37| and Eq. ( |26| ) here). Theoretical error on R(dj/s"f) is 



estimated to be at most a few percent in [37], comparable to the one on ATZ. In hadronic collisions, 
the ratio ATZ is more likely to be measured than R^d'y/sj). 

We also mention here the exclusive radiative decays B — > (p, a;) 7 and B — ► K*^, whose ratios of the 
branching ratios can also be used to determine |T4<i/^ts| UK)- The expected theoretical accuracy on the 
ratio B(B ± — > p ± +~/)/B(B ± —* + 7) is, however, not anticipated to be better than O(20%) Q. 
The corresponding LD-corrections in the ratios of neutral -B-decays, B(B° — > (p°,cj) + ^y)/B(B° — ► 
K*° + 7) are expected to be smaller [ 39 , 40 [ due to their being both colour and (electric)-charge 
suppressed, hence reducing the theoretical uncertainty, but probably not better than ±10%. Finally, 
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Figure 4: The ratio A1Z defined in Eq. (]24), calculated in the HQET-approach, as a function of p 
for three values of r\; solid curve (rj = 0.383), dotted curve (rj = 0.5), dashed curve (77 = 0.27). 



we also note the constraints on |Vfd/Vt s |, which can be obtained from the measurements of the ratios 
of some exclusive two-body non-leptonic decays, such as B(B° — > K*K°)/B(B° — ► (f>K°), advocated in 
Ref. This method may provide interesting results on the CKM ratio, but once data are available 
on the FCNC radiative and semileptonic decays discussed above, they are expected to provide more 
reliable information on the CKM matrix elements Vtd and Vt s - In particular, the ratio A1Z may provide 
one of the most precise determinations of | Vtrf/Vta|. 

We hope that the results presented here will help focus attention on experimental measurements of 
the branching ratios and CP asymmetries in the FCNC decays B — ► (Xj, X s )£ + £~. We also underline 
the need to calculate the next-to-leading order corrections in the CP asymmetries to tame the scale 
dependence. 
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